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A note on mixture representations for the Linnik and Mittag-Leffler 

distributions and their applications* * 

V. Yu. Korolev} A. I. ZeifmarC 


Abstarct: We present some product representations for random variables with the Linnik, Mittag-Leffler 
and Weibull distributions and establish the relationship between the mixing distributions in these representations. 
The main result is the representation of the Linnik distribution as a normal scale mixture with the Mittag-Leffler 
mixing distribution. As a corollary, we obtain the known representation of the Linnik distribution as a scale 
mixture of Laplace distributions. Another corollary of the main representation is the theorem establishing that 
the distributions of random sums of independent identically distributed random variables with finite variances 
converge to the Linnik distribution under an appropriate normalization if and only if the distribution of the 
random number of summands under the same normalization converges to the Mittag-Leffler distribution. 

1 Introduction 

Usually the Mittag-Leffler and Linnik distributions are mentioned in the literature together as examples of 
geometric stable distributions. Since these distributions are very often pointed at as weak limits for geometric 
random sums, there might have emerged a prejudice that the scheme of geometric summation is the only 
asymptotic setting within which these distributions can be limiting for sums of independent and identically 
distributed random variables. This prejudice is accompanied by the suspicion that non-trivial (<5 < 1, a < 2) 
Mittag-Leffler and Linnik laws can be limiting only for sums in which the summands have infinite variances. 

Another obvious reason for which the Mittag-Leffler and Linnik distributions are often brought together in 
the literature is the formal similarity of the Laplace transform of the former and the Fourier-Stieltjes transform 
of the latter. 

The main aim of the present paper is to study the analytic and asymptotic relations between these two 
laws. We will show that actually the link between these two laws is much more interesting. Namely, it turns 
out that the Linnik distribution with parameter a is a scale mixture of the normal distributions with the 
mixing distribution being the Mittag-Leffler law with parameter S = a/2. As a corollary of this result we obtain 
a theorem establishing that the distributions of random sums of independent identically distributed random 
variables with finite variances converge to the Linnik distribution under an appropriate normalization if and 
only if the distribution of the random number of summands under the same normalization converges to the 
Mittag-Leffler distribution. 

Our main tools are mixture representations for the Linnik, Mittag-Leffler and Weibull distributions. Mixture 
representations for the Linnik and Mittag-Leffler laws were the objects of investigation in [3, 5, 6, 16, 23, 19, 20]. 
Some of the results of these papers will be used in what follows and will be presented as lemmas below. As well, 
the proofs of our results are based on some new mixture representations for the Weibull distributions. 

However, we are unaware of the results concerning the possibility of representation of the Linnik distribution 
as a scale mixture of normals. Perhaps, the paper [16] is the closest to this conclusion and exposes the 
representability of the Linnik law as a scale mixture of Laplace distributions with the mixing distribution 
written out explicitly. 

We develop the results of [16] and obtain our main result which is the representation of the Linnik distribution 
as a normal scale mixture with the Mittag-Leffler mixing distribution, thus finding a tight and clear analytical 
link between the Linnik and Mittag-Leffler distributions. As a corollary, we have the representation of the Linnik 
distribution as a scale mixture of Laplace distributions obtained in [16]. Our proof of this result together with the 
explicit formula for the mixing density obtained in [16] lead to a by-product corollary of this representation which 
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is the explicit formula for the distribution density of the ratio of two independent positive strictly stable random 
variables. Another consequence of the main representation is the theorem establishing that the distributions 
of random sums of independent identically distributed random variables with finite variances converge to the 
Linnik distribution under an appropriate normalization if and only if the distribution of the random number 
of summands under the same normalization converges to the Mittag-Lefher distribution. On the one hand, 
this theorem offers a new «asymptotic» link between the Linnik and Mittag-Leffler distributions and, on the 
other hand, it shows that the Linnik law can be limiting for the distribution of sums of independent identically 
distributed random variables with finite variances thus dispelling the suspicion mentioned above. 

The paper is organized as follows. Section 2 presents the definitions and basic properties of the Linnik and 
Mittag-Leffler distributions. Section 3 contains all the definitions and auxiliary results. The proofs of our main 
results are purposely indirect and essentially rely on some new mixture properties of the Weibull distribution 
also presented in Section 3. In Section 4 we prove the representability of the Linnik distribution as the scale 
mixture of normal laws with the Mittag-Leffler mixing distribution and as the scale mixture of the Laplace 
laws with the mixing distribution being that of the ratio of two independent random variables with the same 
strictly stable distribution concentrated on the nonnegative halfline. We use this representation together with 
the result of [16] to obtain a by-product corollary which is the explicit representation of the distribution density 
of the ratio of two independent positive strictly stable random variables. In Section 5 we prove and discuss 
the theorem establishing that the distributions of random sums of independent identically distributed random 
variables with finite variances converge to the Linnik distribution under an appropriate normalization if and 
only if the distribution of the random number of summands under the same normalization converges to the 
Mittag-Leffler distribution. 


2 The Mittag-Leffler and Linnik distributions 

2.1 The Mittag-Leffler distributions 

The Mittag-Leffler probability distribution is the distribution of a nonnegative random variable Mg whose 
Laplace transform is 

Ms) = Ee“ sM * = a £ 0, (1) 

where A>0,0<<5^1. For simplicity, in what follows we will consider the standard scale case and assume 
that A = 1. 

The origin of the term Mittag-Leffler distribution is due to that the probability density corresponding to 
Laplace transform (1) has the form 


fs ML (x) 


1 

x x ~ s 



( -l) n x sn 

o T(Sn + 1) 


-—Eg(-x s ), x ^ 0, 


( 2 ) 


where Eg(z) is the Mittag-Leffler function with index S that is defined as the power series 


w = E 


n=o T(<5n + 1) ’ 


6 > o, zeZ. 


Here T(s) is Euler’s gamma-function, 


TOO 


_.s—1 2 


dz , 


s > 0. 


The distribution function corresponding to density (2) will be denoted Fff' Uj {x). 

With <5 = 1, the Mittag-Leffler distribution turns into the standard exponential distribution, that is, 
F 1 Mi (x) = [1 — e~ x ]l(x ^ 0), x £ R (here and in what follows the symbol 1(C) denotes the indicator function 
of a set C). But with <5 < 1 the Mittag-Leffler distribution density has the heavy power-type tail: from the 
well-known asymptotic properties of the Mittag-Leffler function it can be deduced that 


fs ML (x) 


sin(57r)T((5 + 1) 
7rx s+1 


as x —> oo, see, e. g., [11]- 
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It is well-known that the Mittag-Leffler distribution is stable with respect to geometric summation (or 
geometrically stable). This means that if X\, X 2 , ■ ■ . are independent random variables and N p is the random 
variable independent of X \, X 2 ,... and having the geometric distribution 

P(N p = n) = p(l -p) n_1 , n= 1,2,..., pe(0,1), (3) 

then for each p £ (0,1) there exists a constant a p > 0 such that a p (Xi + ... + X/v p ) => Ms as p —> 0, see, e. 
g., [2] or [13] (the symbol ==>■ hereinafter denotes convergence in distribution). Moreover, as far ago as in 1965 
it was shown by I. Kovalenko [18] that the distributions with Laplace transforms (1) are the only possible limit 
laws for the distributions of appropriately normalized geometric sums of the form a p +... + Xjv p ) as p —>■ 0, 
where X \, X 2 ,... are independent identically distributed nonnegative random variables and N p is the random 
variable with geometric distribution (3) independent of the sequence X\, X 2 , ■ .. for each p £ (0,1). The proofs 
of this result were reproduced in [7, 8] and [9]. In these books the class of distributions with Laplace transforms 
(1) was not identified as the class of Mittag-Leffler distributions but was called class 1C after I. Kovalenko. 

Twenty five years later this limit property of the Mittag-Leffler distributions was re-discovered by A. Pillai 
in [26, 27] who proposed the term Mittag-Leffler distribution for the distribution with Laplace transform (1). 
Perhaps, since the works [18, 7, 8] were not easily available to probabilists, the term class 1C distribution did 
not take roots in the literature whereas the term Mittag-Leffler distribution became conventional. 

The Mittag-Leffler distributions are of serious theoretical interest in the problems related to thinned (or 
rarefied) homogeneous flows of events such as renewal processes or anomalous diffusion or relaxation phenomena, 
see [29, 10] and the references therein. 

2.2 The Linnik distributions 

In 1953 Yu. V. Linnik [22] introduced the class of symmetric probability distributions defined by the characteristic 
functions 

V>a(t) = t £ R, (4) 

where a £ (0,2]. Later the distributions of this class were called Linnik distributions [17] or a-Laplace distributions 
[25]. In this paper we will keep to the first term that has become conventional. With a = 2, the Linnik distribution 
becomes the Laplace distribution corresponding to the density 

f A {x) = i e -H, ieR. (5) 

A random variable with Laplace density (5) and its distribution function will be denoted A and F A (x), 
respectively. 

The Linnik distributions possess many interesting analytic properties such as unimodality [21] and infinite 
divisibility [3], existence of an infinite peak of the density for a ^ 1 [3], etc. However, perhaps, most often Linnik 
distributions are recalled as examples of geometric stable distributions. 

A random variable with the Linnik distribution with parameter a will be denoted L a . Its distribution 
function and density will be denoted F A and f A , respectively. As this is so, from (4) and (5) it follows that 
F£{x) = F A {x), i£l. 


3 Basic notation and auxiliary results 


Most results presented below actually concern special mixture representations for probability distributions. 
However, without any loss of generality, for the sake of visuality and compactness of formulations and proofs 
we will formulate the results in terms of the corresponding random variables assuming that all the random 
variables mentioned in what follows are defined on the same probability space (f2, 21, P). 

The random variable with the standard normal distribution function $(x) will be denoted X , 


P(X < x) = <fr(x) 



x £ R. 


Let d^rr), x £ R, be the distribution function of the maximum of the standard Wiener process on the unit 
interval, T(x) = 2$(max{0,^}) — 1, x £ R. It is easy to see that \k(x) = P(|X| < x). Therefore, sometimes 
T(a;) is said to determine the half-normal or folded normal distribution. 

Throughout the paper the symbol = will denote the coincidence of distributions. 
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The distribution function and density of the strictly stable distribution with the characteristic exponent a 
and shape parameter 6 defined by the characteristic function 

fa,e(t) = exp { — |f| a exp{—iRrflasignf}}, let, ( 6 ) 

with 0 < a < 2, \6\ ^ min{l, ^ — 1 }, will be denoted by G a ^(x) and g a ,e(x), respectively (see, e. g., [30]). Any 
random variable with the distribution function G ay g(x) will be denoted S a , 9 - 

From ( 6 ) it follows that the characteristic function of a symmetric (9 = 0) strictly stable distribution has 
the form 

fa,o (i) = e- |i|0 , fei. (7) 

From (7) it is easy to see that S 2 ,o = \[2X. 

Lemma 1. Let a £ (0,2], a' £ (0,1]. Then 


S, 


otot r ,0 


jf „ nl/a 


where the random variables on the right-hand side are independent. 

PROOF. See, e. g., theorem 3.3.1 in [30]. 

Corollary 1 . A symmetric strictly stable distribution with the characteristic exponent a is a scale mixture 
of normal laws in which the mixing distribution is the one-sided strictly stable law (9 = 1) with the characteristic 
exponent a/ 2 : 

nOO 

G a , o( x )= $(x/V2z)dG a / 2 ,i{z), xeR. 

Jo 

In terms of random variables the statement of corollary 1 can be written as 


5q:,0 


= X^2S a/2A 


( 8 ) 


with the random variables on the right-hand side being independent. 

Let 7 > 0. The distribution of the random variable W 7 : 

P(W 7 <ar) = [l-e“ :c 7 ]l(x>0), iel, 

is called the Weibull distribution with shape parameter 7 . It is obvious that W\ is the random variable with the 
standard exponential distribution: P(W\ < x) = [l — e -x ] l(x ^ 0). The Weibull distribution with 7 = 2, that 
is, P(W 2 < x) = [l — e~ x2 \ 1 (x ^ 0) is called the Rayleigh distribution. 

It is easy to see that if 7 > 0 and 7 ' > 0, then P(W ^/ 7 ^ x) = P(W 7 < ^ x J ) = e ~ x77 = P(W 77 / ^ x), 
x ^ 0 , that is, for any 7 > 0 and 7 ' > 0 

W 7 y = W 7 V 7 . (9) 

In particular, for any 7 > 0 we have W 1 = . 

It can be shown that each Weibull distribution with parameter 7 £ (0,1] is a mixed exponential distribution. 
In order to prove this we first make sure that each Weibull distribution with parameter 7 £ (0,2] is a scale 
mixture of the Rayleigh distributions. 

For a £ ( 0 , 1 ] denote V a = S~\, where S a ,i is a random variable with one-sided strictly stable density 
9a, i(x). 

Lemma 2. For any 7 £ (0, 2] we have 

w 7 ± W 2y Jv^~ 2 , (10) 

where the random variables on the right-hand side of (10) are independent. 

PROOF. Write relation (8) with a replaced by 7 in terms of characteristic functions with the account of (7): 

poo 

e -|t | 7 = / exp{— t 2 z}g 7 / 2 ,i(z)dz, t £ R. ( 11 ) 

J 0 

Formally letting |t| = x in (11), where x ^ 0 is an arbitrary nonnegative number, we obtain 

pOO 

P(W 7 > x) = e -x7 = / exp{— x 2 z}g 1 / 2 ,i(z)dz. ( 12 ) 

Jo 
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( 13 ) 


At the same time it is obvious that if W 2 and SL/ 2,1 are independent, then 

P{W 2 ^V l/2 > x) = P (W 2 > XyJS 7 / 2 ,i) = J exp{-x 2 z}g l/ 2 ,i(z)dz. 

Since the right-hand sides of (12) and (13) coincide identically in x ^ 0, the left-hand sides of these relations 
coincide as well. The lemma is proved. 

Lemma 3. For any 7 £ (0,1], the Weibull distribution with parameter 7 is a mixed exponential distribution: 


W 1 = W^j. (14) 

where the random variables on the right-hand side of (14) are independent. 

PROOF. It is easy to see that P{Wfz x) = P(W\ ^ x 1 ) = e ~ x ’ 7 = P(W 7 > x), x ^ 0, that is, 

W 1 = W\ h . (15) 

From (15) it follows that W 2 = \JW\. Therefore, from lemma 2 it follows that for 7 £ (0, 2] we have 

w 7 = W 2 Jv ^2 = yJw ^/2 

or, with the account of (15), 

W 7/2 £ W 2 = WlVr^/2- 

Re-denoting 7/2 1 —> 7 £ (0,1], we obtain the desired assertion. 

In [3] the following statement was proved. Here its formulation is extended with the account of (9). 

Lemma 4 [3]. For any a £ (0, 2], the Linnik distribution with parameter a is a scale mixture of a symmetric 
stable distribution with the Weibull mixing distribution with parameter a/2, that is, 

L a = S afi W a = S afi Jw^2, 


where the random variables on the right-hand side are independent. 

Lemma 5. For any 5 £ (0,1], the Mittag-Leffler distribution with parameter S is a scale mixture of a 
one-sided stable distribution with the Weibull mixing distribution with parameter 8/2, that is, 


Ms = Ss,lW B = Ss,iy/Wsp t 

where the random variables on the right-hand side are independent. 

PROOF. This statement has already become folklore. For the purpose of convenience we give its elementary 
proof without any claims for priority. Let Ss.i be a positive strictly stable random variable. As is known, its 
Laplace transform is = Ee -55,5 ’ 1 = e~ s , s ^ 0. Then with the account of (9) by the Fubini theorem the 
Laplace transform of the product Sy 1 Ws is 


pOO 

Eexp{—sSyiWa} = Eexpl-s^iITi 7,5 } = EE(exp{-s5 a ,iIF 1 1 / 5 }|lTi) = / e“ ( szV ^ e~ z dz = 

Jo 


e~ z{sS+l) dz = r = Ee ~ sMs , s > 0. 


jo 1 + s 

The lemma is proved. 

Let p £ (0,1). In [19] it was demonstrated that the function 

sin(7rp) 


/,(*) = 


x £ ( 0 , 00 ), 


Ttp[x 2 + 2xcos{j:p) + 1 ] 

is a probability density on (0, 00 ). Let K p be a random variable with density (16). 
Lemma 6 [19]. Let 0 < 8 < 8' ^ 1 and p = 8/8' < 1. Then 


( 16 ) 


Ms = Ms'K l p /s 
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where the random variables on the right-hand side are independent. 

With S' = 1 we have 

Corollary 2 [19]. Let 0 < S < 1. Then the Mittag-Leffler distribution with parameter S is mixed exponential. 

M s = K] ,S W 1 


where the random variables on the right-hand side are independent. 
Let 0 < a < a' ^ 2. In [16] it was shown that the function 


•WO) = 


a' smfna/a')x" 




7r[l + x 2a + 2x a cos( 7 ra/a')] ’ 


x > 0 , 


(17) 


is a probability density on ( 0 , oo). Let Q a ,a' be a random variable whose probability density is a '( x )- 
Lemma 7 [16]. Let 0 < a < a' ^ 2. Then 


L n = L n 


Qo 


where the random variables on the right-hand side are independent. 

With a! = 2 we have 

COROLLARY 3 [16]. Let 0 < a < 2. Then the Linnik distribution with parameter a is a scale mixture of 
Laplace distributions corresponding to density (5): 

L a = AQ q , 2 


where the random variables on the right-hand side are independent. 

For the sake of completeness, we will demonstrate that the Weibull distributions possess the same property 
as the Linnik and Mittag-Leffler distributions presented in lemmas 6 and 7: any distribution of the corresponding 
class can be represented as a scale mixture of a distribution from the same class with larger parameter. 

Relation (14) implies the following statement generalizing lemmas 2 and 3 and stating that the Weibull 
distribution with an arbitrary positive shape parameter 7 is a scale mixture of a Weibull distribution with an 
arbitrary positive shape parameter 7 ' > 7 . 

Lemma 8 . Let 7 ' > 7 > 0 be arbitrary numbers. Then 

W 7 = Wy ' Vl h \ 


where a = 7 /y' £ ( 0 , 1 ) and the random variables on the right-hand side are independent. 

PROOF. In lemma 3 we showed that a Weibull distribution with parameter a £ (0,1] is a mixed exponential 
distribution. Indeed, from (14) it follows that 

pQO 

e~ x = P (W a > x) = P(Wi > S a ,\x) = / e~ zx g a ,i(z)dz, x ^ 0. 

Jo 

Therefore, for any 7 ' > 7 > 0, denoting a = 7 / 7 ' (as this is so, a £ (0, 1)), for any i£lwe obtain 


P(W 7 > x) = e~ xl = e- xl ' a = P (W a > x 7 ') = P(Wi > S a , ix 7 ') = 


pOO , pOO 

/ e~ zxl g a ,i{z)dz = / P(Wy > xz lh )g a ,i(z)dz = P(Wy • I4 1/7 > x), 
Jo Jo 


The lemma is proved. 

It should be noted that if 0 < 7 < 7 ' < 2, then the assertion of lemma 8 directly follows from theorem 3.3.1 
of [30] due to the formal coincidence of the characteristic function of a strictly stable law and the complementary 
Weibull distribution function (see the proof of lemma 2). 
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4 Representation of the Linnik distribution as a scale mixture of 
normal or Laplace distributions and related results 

In all the products of random variables mentioned below the multipliers are assumed independent. 

The following statement presents the main result of this paper. 

Theorem 1. Let a £ (0,2], a' £ (0,1]. Then 

L aa . = S Oi0 M^ Q 


PROOF. From lemma 4 we have 


— Q 

— 0 



Continuing (18) with the account of lemma 1, we obtain 


L 


aa' 


d_ c qI/c* 



From (9) and lemma 5 it follows that 


C 1 /" 

‘-V,i 



— Q l / a W / 
** acx ' 


= (Sa^lW^ = M^°. 


(18) 


(19) 


The theorem is proved. 

As far as we know, the following result has never been explicitly presented in the literature in full detail 
although the property of the Linnik distribution to be a normal scale mixture is something almost obvious. 

Corollary 4. For each a £ (0, 2], the Linnik distribution with parameter a is the scale mixture of zero-mean 
normal laws with mixing Mittag-Leffler distribution with twice less parameter a/2: 


L a = X^2M a/2 , 

where the random variables on the right-hand side are independent. 

It should be noted that in this case representation (19) takes the form 


L a — X ^j2S a /2sW a /2- 


( 20 ) 


From lemma 3 we have 


Hence, from (20) it follows that 



Wi 


L a =X 



S„ 


/2,1 

/2,1 


where the independent random variables S a / 2 p and S' a ^ 2 1 have one and the same one-sided strictly stable 
distribution with characteristic exponent a/2 and are independent of the exponentially distributed random 
variable W\. It is well known that 


Xsfrwl = A 


( 21 ) 


(see, e. g., the example on p. 272 of [1]). Therefore we obtain one more mixture representation for the Linnik 
distribution. 

Theorem 2. For each a £ (0,2], the Linnik distribution with parameter a is the scale mixture of the Laplace 
laws corresponding to density (5) with mixing distribution being that of the ratio of two independent random 
variables having one and the same one-sided strictly stable distribution with characteristic exponent a/2: 


L„ = A, 


IS, 


a/2,1 


S" 


a/2,1 


where the random variables on the right-hand side are independent. 
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It is easy to see that scale mixtures of Laplace distribution (5) are identifiable, that is, if 


A Y = AY' 

where Y and Y are nonnegative random variables independent of A, then Y = Y'. Indeed, with the account of 
(21), the last relation turns into 

X^j2W x Y 2 = XsjTW x ifY') 2 , (22) 

where the random the multipliers on both sides are independent. But, as is known, scale mixtures of zero-mean 
normals are identifiable (see [28]). Therefore, (22) implies that 

W{Y 2 = Wi(F') 2 - (23) 


The complementary mixed exponential distribution functions of the random variables related by (23) are the 
Laplace transforms of Y 2 and (Y 1 ) 2 , respectively. Relation (23) means that these Laplace transforms identically 
coincide: 


z dP(Y‘ 


nOO 

<z)= e~ sz dP{(Y') 2 < z), 
Jo 


s^O. 


Hence, the distributions of the random variables Y 2 and ( Y ') 2 coincide and hence, the distributions of Y and 
Y' coincide as well since these random variables were originally assumed nonnegative. 

Comparing the statement of theorem 2 with the assertion of corollary 3 with the account of identifiability 
of scale mixtures of Laplace distributions (5) we arrive at the relation 


Qa,2 


d 



(24) 


The combination of (17) and (24) gives one more, possibly simpler, proof of the following by-product result 
concerning the properties of stable distributions obtained in [4]. This result offers an explicit representation 
for the density of the ratio of two independent stable random variables in terms of elementary functions 
although with the exception of one case, the Levy distribution (a = 1), such representations for the densities 
of nonnegative stable random variables themselves do not exist. 

COROLLARY 5. Let S a ,i and, S' a 1 be two independent random variables having one and the same one-sided 

strictly stable distribution with characteristic exponent a £ (0,1). Then S a ,\/S' a x = Kf/ a = Q 2 a 2 , that is, the 
probability density p a (x) of the ratio S at i/S' a x has the form 


Pa (x) = f® i(x) 


sin(7ra)a:“ 1 

7r[l + x 2a + 2x a cos(7ra)] ’ 


x > 0. 


The proof immediately follows from the observation that p a {x) = (2y/x) 1 f^a 2 (V%) = fa i( x )- 
As concerns the Mittag-Leffler distribution, from lemmas 3 and 5 we obtain the following statement 
analogous to theorem 2. 

Theorem 3. For each S € (0,1], the Mittag-Leffler distribution with parameter 5 is the mixed exponential 
distribution with mixing distribution being that of the ratio of two independent random variables having one and 
the same one-sided strictly stable distribution with characteristic exponent 5: 

*6,1 


where the random variables on the right-hand side are independent. 

From theorem 3 and corollary 5 we obtain the following representation of the Mittag-Leffler distribution 
function F^ L (x): 


? ML 


(x) = 1 - 


sin(7r<5) 


z 6-i e -zx dz 


1 + z 2S + 2 z s cos(ttS) ’ 


x > 0. 


(25) 


Representation for the Linnik distribution similar to (25) was obtained in [19]. 

Using lemmas 1, 4 and 5 it it is possible to obtain more product representations for the Mittag-Leffler- and 
Linnik-distributed random variables and hence, more mixture representations for these distributions. 









5 Convergence of the distributions of random sums to the Linnik 
distribution 

Product representations for the random variables with the Linnik and Mittag-Leffler distributions obtained 
in the previous works were aimed at the construction of convenient algorithms for the computer generation 
of pseudo-random variables with these distributions. The mixture representation for the Linnik distribution 
as a scale mixture of normals obtained in corollary 4 opens the way for the construction in this section of a 
random-sum central limit theorem with the Linnik distribution as the limit law. Moreover, in this version of 
the random-sum central limit theorem the Mittag-Leffler distribution must be the limit law for the normalized 
number of summands. 

Recall that the symbol =>• denotes the convergence in distribution. 

Consider a sequence of independent identically distributed random variables X\ , X^, ■ . . defined on the 
probability space (f 1, 21, P). Assume that EXi = 0, 0 < a 2 = DX± < oo. For n £ N denote S* = X\ + ... + X n . 
Let Zi, Z 2 , ■ ■ ■ be a sequence of integer-valued nonnegative random variables defined on the same probability 
space so that for each n ^ 1 the random variable Z n is independent of the sequence X±, X 2 , ■ ■ . For definiteness, 
in what follows we assume that ]Cj=i = 0- 

p 

Recall that a random sequence Z\, Z2, ■ • ■ is said to infinitely increase in probability {Z* n —y 00), if P (Z n ^ 
to) —y 0 as n —y 00 for any to £ ( 0 , 00). 

The proof of the main result of this section is based on the following version of the random-sum central 
limit theorem. 

Lemma 10 . Assume that the random variables Xi, X2, ■ ■ ■ and Z\,Z2,-.. satisfy the conditions specified 

p 

above and, moreover, let Z n —> 00 as n —y 00. A distribution function F(x) such that 



as n —y 00 exists if and only if there exists a distribution function H ( x) satisfying the conditions 

H( 0) = 0, F{x) = *(-^)dff(y), ieR, 


and P(Z n < nx) =>• H{x) (n —► 00 ). 

PROOF. This statement is a particular case of a result proved in [14], also see theorem 3.3.2 in [9]. 

The following theorem gives a criterion (that is, necessary and sufficient conditions) of the convergence of 
the distributions of random sums of independent identically distributed random variables with finite variances 
to the Linnik distribution. 

Theorem 4. Let a £ (0,2]. Assume that the random variables X\,X 2 ,. ■ ■ and Z\,Z 2 ,... satisfy the 

p 

conditions specified above and, moreover, let Z n — > 00 as n — > 00 . Then the distributions of the normalized 
random sums S* Zn converge to the Linnik law with parameter a, that is, 

Vayn / 

as n —> 00 , if and only if 

— => 2¥„/ 2 (n —y 00 ). 
n 

PROOF. This statement is a direct consequence of corollary 4 and lemma 10 with H(x) = F^j^(x/2). 

The convergence of the distributions of the normalized indices to the Mittag-Leffler distribution Fg IL 

is the main condition in theorem 4. Now we will give two examples of the situation where this condition can 
hold. The first example is trivial and is based on the geometric stability of the Mittag-Leffler distribution. The 
second example relies on a useful general construction of nonnegative integer-valued random variables which, 
under an appropriate normalization, converge to a given nonnegative (not necessarily discrete) random variable, 
whatever the latter is. 
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Example 1. Let <5 E (0,1) be arbitrary. For every n € N let N\/ n be a random variable having the geometric 
distribution (3) with P = ~ independent of the sequence Vi,^,... of independent identically distributed 
nonnegative random variables such that 


n 


- 1/5 


2 -^ 3 =1 


Yi 


2 M s 


(26) 


asn-> oo. To provide (26), the distributions of the random variables Yi, I 2 , ■ ■ • should belong to the domain of 
the normal attraction of the one-sided strictly stable law with characteristic exponent S. As Z n for each n € N 
take 


Z„. — 


,1-1 /sy N ^ Y 
^3= 1 


where square brackets denote the integer part. Then 


^ = n- 1 ' 5 
n 


Z^ J= i 


Yj 



(27) 


where curly braces denote the fractional part. Since the second term on the right-hand side of (27) obviously 
tends to zero in probability, from (26) it follows that Z n jn ==> 2 Ms asn-> 00. 

Example 2. In the book [9] it was proposed to model the evolution of non-homogeneous chaotic stochastic 
processes, in particular, the dynamics of financial markets by compound doubly stochastic Poisson processes 
(compound Cox processes). This approach got further grounds and development, say, in [1, 15]. According to 
this approach the flow of informative events, each of which generates the next observation, is described by the 
stochastic point process P(Y(t)) where P(t), t > 0, is a homogeneous Poisson process with unit intensity and 
Y(t), t > 0, is a random process independent of P(t) and possessing the properties: Y (0) = 0, P(Y (t) < 00) = 1 
for any t > 0, the trajectories Y(t) are non-decreasing and right-continuous. The process P(Y(t)), t > 0, is 
called a doubly stochastic Poisson process (Cox process) [12]. 

Within this model, for each t the distribution of the random variable P{Y(t)) is mixed Poisson. For vividness, 
consider the case where in the model under consideration the parameter t is discrete: Y(t) = Y(n) =Y n ,n&N, 
where {Y n } n ^ 1 is an infinitely increasing sequence of nonnegative random variables such that Y n+ i(co) ^ Y n {ui) 
for any u € 9, rt ^ 1. Here the asymptotics n —> 00 may be interpreted as that the intensity of the flow of 
informative events is (infinitely) large. 

Assume that the random variable Ms is independent of the standard Poisson process P(t), t ^ 0. For each 
natural n take Y n = 2 nM$. Respectively, let Z n = P{Y n ) = P{2nMg), n ^ 1. It is obvious that the random 
variable Z n so defined has the mixed Poisson distribution 


P(Z n = k) = P(P( 2 nMs) = k) = j e- 2 nz ^^dF s ML {z) k = 0,1,... 

This random variable Z n can be interpreted as the number of events registered up to time n in the Poisson 
process with the stochastic intensity distributed as 2 Ms. 

Denote A n (z) = P (Z n < 2 nz), z ^ 0 (A n (z) = 0 for 2 < 0 ). It is easy to see that A n (z) =>• F^ L {z) as 
n —> 00. Indeed, as is known, if n(a;; £) is the Poisson distribution function with the parameter £ > 0 and E(x\ c) 
is the distribution function with a single unit jump at the point c E R, then H(£x',£) => E{x\ 1 ) as l —> 00. 
Since for i£l 

pOO 

A n (x) = / H( 2 nx; 2 nz)dFg IL (z), 

Jo 

then by the Lebesgue dominated convergence theorem, as n —> 00, we have 

poo px 

A n (x) =► / E(x/z-, 1 )dF s ML (z) = / dF s ML (z) = F s ml (x), 

Jo Jo 

p 

that is, the random variables Z n defined above satisfy the condition of theorem 4 . Moreover, Z n —> 00 as 
n —> 00 since P (Ms = 0 ) = 0 . 

The authors thank Dmitri Alekritskii who carefully read the draft text of the paper and pointed out some 
inaccuracies. 
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